Anderson transition in a three dimensional kicked rotor 
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We investigate Anderson localization in a three dimensional (3d) kicked rotor. By a finite size 
scaling analysis we have identified a mobility edge for a certain value of the kicking strength k = k c . 
For k > k c dynamical localization does not occur, all eigenstates are delocalized and the spectral 
correlations are well described by Wigner-Dyson statistics. This can be understood by mapping the 
kicked rotor problem onto a 3d Anderson model (AM) where a band of metallic states exists for 
sufficiently weak disorder. Around the critical region k « k c we have carried out a detailed study of 
the level statistics and quantum diffusion. In agreement with the predictions of the one parameter 
scaling theory (OPT) and with previous numerical simulations of a 3d AM at the transition, the 
number variance is linear, level repulsion is still observed and quantum diffusion is anomalous with 
(p 2 ) oc t 2//3 . We note that in the 3d kicked rotor the dynamics is not random but deterministic. 
In order to estimate the differences between these two situations we have studied a 3d kicked rotor 
in which the kinetic term of the associated evolution matrix is random. A detailed numerical 
comparison shows that the differences between the two cases are relatively small. However in the 
deterministic case only a small set of irrational periods was used. A qualitative analysis of a much 
larger set suggests that the deviations between the random and the deterministic kicked rotor can 
be important for certain choices of periods. Contrary to intuition correlations in the deterministic 
case can either suppress or enhance Anderson localization effects. 
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The quantum kicked rotor [l], Q , namely, a particle in 
a circle periodically kicked (with period T) by a smooth 
potential V(q) = kcos(q), 

H = ^ + V(q)J2S(t-Tn), (1) 

n 

has played a central role in the development of quan- 
tum chaos. This is hardly surprising since, despite its 
simplicity, it has a very rich and highly non-trivial dy- 
namics. Moreover, it can be studied experimentally [1, 0| 
by using cold atoms techniques. 

The classical dynamics undergoes a gradual crossover 
from quasi-integrable to fully chaotic as the kicking 
strength k is increased. Therefore it is possible to study 
in detail all types of intermediate behavior between chaos 
and integrability. Moreover it is specially suited for nu- 
merical simulations since the evolution matrix in a basis 
of plane waves has a simple form. This was important 
in the early days of quantum chaos when numerical sim- 
ulations in conservative 2d chaotic systems were quite 
expensive. 

The quantum dynamics of the kicked rotor is also of 
interest. For short time scales and k sufficiently large 
quantum and classical motion are diffusive. However for 
longer time scales numerical simulations show [l| that 
quantum diffusion in momentum space eventually stops, 
namely, limt_ >oc (p 2 (i)) tends to a constant. By contrast 
classically (p 2 (t)) oc t for all t. This rather counterintu- 



itive feature, usually referred to as dynamical localization 
, was fully understood || after mapping the kicked ro- 
tor problem onto a Id AM with a pseudo-random poten- 
tial where localization is well established. In our model 
the role of the mean free path is played by the kicking 
strength k. The mapping onto a Id AM is carried out in 
two steps: first one writes down the evolution operator 
U associated to Eq.JTj, 

-ip 2 T iV(q) 

U — e 2h e h . (2) 
In a basis of plane waves, U has matrix elements, 

U lm = (-i) l -™ e - im2 ^ Jj_ m (3) 

where J n (k) is a Bessel function of first kind. Then a 
given eigenstate of U with eigenvalue exp(— iu) can be 
transformed into an effective Hamiltonian problem [||, 

7-Lip n = e„-0n + F(m - n)tp m (4) 

m 

where e„ = tan(w/2 — rn 2 /4) and 

F(m-n) = -— / dqtan(V(q)/2)e-^ m ^ 

2n J — 7T 

oc exp(— A\m — n\) (5) 

with A a decreasing function of k and r = T/h. For 
comparison we note that the standard Id AM heavily 
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used in numerical and analytical studies of localization is 
given by, 

T-tlpn = e n 1p n + 1p n+ i + lp n -i (6) 

where e n are independent random numbers extracted 
from a Gaussian or box distribution. The quantum dy- 
namics of the Hamiltonian Eq.((6|) is well understood 
@) S eigenstates are exponentially localized and 
quantum diffusion eventually is arrested for any amount 
of disorder. 

We note that strictly speaking there are important 
differences between Eq.© and Eq.(|3}. For instance in 
Eq.((4|) the diagonal disorder e n is deterministic and the 
hopping is short range but it is not restricted to nearest 
neighbors as in Eq.©. The latter is not really relevant 
as it is known [9j that exponential localization for all en- 
ergies and disorder persists provided that the hopping 
term decays faster than l/|n — m\. Potentials in Eq.JI]) 
leading to a decay of F(m — n) ~ l/\n — m\ in Eq.Q 
induce a metal-insulator transition even in Id ddif. 

The non-randomness of the diagonal disorder in Eq. ^ 
is, at least superficially, a more serious deviations from 
Eq.©. However, in Ref. [TTI] it was showed that provided 
t ^> 1 is irrational with a sufficiently poor approxima- 
tion by a rational number, the potential was ergodic and 
to a great extent not very different from a true random 
potential. This seems to be related to the fact that in the 
n — > oo limit all derivatives of e„ are not bounded. For 
the sake of completeness we mention that a) it was not 
until recently [12J that a rigorous proof of localization for 
the original model Eq.([T]), not involving a mapping onto 
an AM, is available; and b) localization in more general 
time dependent potentials has been recently studied in 
the mathematical literature. The main conclusion [l3j 
is that if the potential is periodic and smooth in time, 
dynamical localization will generally occur. 

A natural question to ask is to what extent these re- 
sults hold in higher dimensions. One of the more strik- 
ing predictions of localization theory 0, [§[ , corroborated 
by numerical simulations, is that eigenstates of the 2d 
version of Eq.® are still exponential localized for any 
disorder and energy. As in the Id case, quantum diffu- 
sion eventually stops as well. However in more than 2d a 
metal insulator transition occurs for a critical value of the 
disorder. The critical disorder increases with the spatial 
dimensionality of the system [1, 0, Q . It is thus highly 
desirable to determine a) whether it is possible to extend 
the above mapping to higher dimensions, and b) whether 
the resulting extension of Eq. (0J still follows the predic- 
tions of localization theory despite the potential being 
non-random. 

We note that if the two models are still similar higher 
dimensional kicked rotors could be used to study exper- 
imentally the Anderson transition by using cold atoms 
techniques. This is specially relevant as experimental 
detection of localization has proved to be quite elusive. 
Moreover it is also important to understand in detail how 
deviations from an ideal random distribution affects the 



quantum dynamics of the system. After all any realistic 
random potential will always present some sort of devia- 
tion from a given ideal random distribution. 

In fact there are already in the literature different 
propositions to map a generalization of Eq.© onto a 
higher dimensional AM. In Ref. [14| a multifrequency, 
time dependent potential 

V(q) = kcos(q + Uit + u 2 t) (7) 

was mapped onto a 3d version of Eq. © ; namely, an AM 
with short range hopping and diagonal disorder 

£n, ni ,n 2 ~ tan(rn 2 + LUiTli + uj 2 n 2 ) ■ (8) 

with n, n\,n2 integers numbers. A numerical analysis 
of the quantum diffusion [14j | showed that, as in a 3d 
AM with a truly random potential, the system avoids 
dynamical localization for k > k c . Recently it has been 
possible to model this generalized kicked rotor by using 
cold atoms in an optical lattice [l5[. In agreement with 
previous results no dynamical localization was observed 
for fc > fc c . 

However we note that in the rii and n 2 direction the 
resulting diagonal disorder is quasi-periodic rather than 
quasi-random. It is thus unclear whether this type of 
kicked rotor can be used as an effective model to study 
a metal-insulator transition induced by disorder. Indeed 
the value of the critical exponent [16[ controlling the 
divergence of the localization length in generalizations of 
Eq. ((SJ) with d > 3 frequencies w\ . . .Wd sharply disagrees 
with numerical results [1, [TtJ of the d-dimensional 
equivalent of Eq. © . 

Another generalized kicked rotor with potential, 

V(qi,q2) = fccos(<7i)cos(<7 2 ) (9) 

was studied in Ref. [I8J . In this case the mapping is onto 
a 2d AM with short range hopping and diagonal disorder, 

£m,na ~ tannin? + r 2 n\). (10) 

with ni,n 2 integers and ti,t 2 irrationals. This diagonal 
pseudo-disorder is a natural extension of the one dimen- 
sional case Eq. © and consequently shares many of the 
properties of a true random potential. For instance, the 
dynamical localization occurs for any value of the kicking 
strength, and the localization length increases exponen- 
tially with fc. 

In this paper we study a 3d generalization of the above 
model, namely, a 3d kicked rotor with 

V{qi,Q2,qa) = k cos(q 1 ) cos(q 2 ) cos(q 3 ) . (II) 

It will be shown that this model maps onto a 3d version of 
Eq.© with short range hoping and an on-site potential 
that shares many quasi-random features of the Id and 2d 
cases. 

By a careful finite size scaling analysis it will shown 
that this model undergoes a metal insulator transition at 
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a certain kicking strength k — k c . For technical reasons 
we have investigated in detail only a small sets of irra- 
tional periods {t;} with i — 1,2,3. Within this set the 
numerical value of k c depends weakly on the {t^}. Dy- 
namical localization consequently only occurs for k < k c . 
We then study in detail level statistics and quantum dif- 
fusion in the critical region k ~ k c . In agreement with 
previous numerical results of the level statistics of the 3d 
AM at the transition, we have found that a) the spectral 
rigidity is linear; b) as in a metal, level repulsion is still 
observed and, c) as in an insulator, the asymptotic decay 
of the level spacing distribution is exponential. Moreover, 
in agreement with the the one parameter scaling theory 
(OPT), the diffusion is anomalous at the transition with 
(p 2 ) oci 2 / 3 . 

Finally we investigate to what extent Anderson local- 
ization is affected by the pseudo random nature of the 
potential. In order to proceed we compare the results 
of previous sections with those coming from a 3d kicked 
rotor in which the phase of the kinetic term of the evo- 
lution matrix is random. The main differences between 
the random and pseudo random cases are as follow: a) 
the critical kicking strength k c is slightly larger (k c ~ 2.4 
versus k c ~ 2.3) in the random case; b) the slope of the 
spectral rigidity is smaller in the random case and c) sta- 
tistical fluctuations are much stronger in the pseudo ran- 
dom case. These results are consistent with the heuristic 
picture that localization will be stronger the more ran- 
dom/uncorrelated the potential is. However a qualita- 
tive analysis of a much larger set of {n} suggests that 
the deviations between the random and the determinis- 
tic kicked rotor can be important for certain choices of 
periods. Contrary to intuition certain choices of {rt} can 
enhance Anderson localization effects and lead to a larger 
k c than in the random case. 

We start with a brief review of the features expected 
at the metal- insulator transition according to numerical 
simulations and OPT 0. 



A. The one parameter scaling theory and 
properties of the 3d metal-insulator transition 

In the context of disordered systems the one parameter 
scaling theory (OPT) 0] provides a valuable framework 
to understand localization effects. In two and lower di- 
mensions it predicts the exponential localization of the 
eigenstates and the arrest of diffusion for any amount 
of disorder [3]. In three and higher dimension a metal 
insulator transition is expected for a critical amount of 
disorder. Thus for disorder below the critical one the 
wavefunctions are extended through the sample. In the 
opposite limit, wave functions are exponentially localized 
as in d < 2. 

A key concept in this theory is the dimensionless con- 
ductance g introduced by Thouless [l!|. It is defined as 
g = Et/A where, Et, the Thouless energy, is an energy 
scale related to the diffusion time to cross the sample and 



A oc l/L d is the mean level spacing. In the diffusive limit 
Et oc 1/L 2 and therefore g oc L d ~ 2 . On the other hand 
if the the particle is exponentially localized, g oc e~ L /£ 
where £ is the localization length and L is the system 
linear size. 

With this simple input it is possible to analyze local- 
ization effects in 3d: For a metal limx,_ ) . 00 g{L) — > oo. For 
an insulator liniL^oo g{L) — > 0. Since g(L) is continuous 
and monotonous, the metal insulator transition is charac- 
terized by a size independent conductance g — g c . Since 
g c does not depend on the system size and A oc 1/L 3 , 
the Thouless energy must scale as Et oc 1/L 3 . This only 
can happen if the diffusion is anomalous at the transition 
with 



(r 2 }~t 



2/3 



(12) 



Anomalous diffusion is one of the main predictions of the 
OPT about the metal insulator transition. 

Additional information about the transition comes 
mainly from numerical calculations [l7j . Level statistics 
[20L l2l| are intermediate between those of a metal and 
those of an insulator. For instance level repulsion is ob- 
served as in a disordered metal. By contrast long range 
spectral correlations such as the spectral rigidity A3 (A) 
are closer to that of an insulator [22j, |23j . For A>1, 



A 3 (A) 



A 4 



(A 3 - 2A 2 ir + x 3 )£ 2 (x)dx - xA/15 



where % < 1 at the metal insulator transition and x = 1 
for an insulator, A is a spectral window containing A 
eigenvalues in units of the mean level spacing A, S 2 (A) = 
(NjD — (N\) 2 is the number variance (JVa is the number 
of eigenvalues in an interval of length A) . 

The critical exponent controlling the divergence of the 
localization length at the transition is supposed to be a 
decreasing function of the space dimensionality. In 3d 
it is approximately given by v « 3/2 17] . Recently one 
of us showed [22j that some of these features can be ob- 
tained analytically by combining the OPT with the self- 
consistent theory of localization 24 1 . 



I. CALCULATION OF THE CRITICAL 
STRENGTH k c AND CRITICAL EXPONENT v 

We study the quantum dynamics of a 3d kicked rotor 
with a smooth potential: 

1 



+V(q 1 ,q 2 ,q 3 )^2S(t - n) 



(13) 



with 



V(qi,q 2 ,q 3 ) = k cos(gi) cos(g 2 ) cos(g 3 ), (14) 



T = h = 1, Ti,T2,T3 are irrational numbers not close to 
any rational multiple of n. It is important to choose {r^} 
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FIG. 1: Scaling variable 77 given by Eq. (|19p as a function of k for different system sizes N. (a) and (b) are the results for two 
different sets of {n} generated with Method 1 (see text), (c) are the averaged results of (a) and (b). The finite size scaling 
analysis suggests that a metal-insulator transition occurs at k c = 2.27 ± 0.07 with a critical exponent v = 1.61 ± 0.25. Plots 
(d), (e) and (f) are the equivalent of plots (a), (b) and (c) but with {n} obtained by Method 2 instead. In this case we obtain 
k c = 2.35 ± 0.1 with v = 1.67 ± 0.27. 



in such a way that correlations caused by the proximity 
to a rational number are suppressed. However it is hard 
to quantitatively assess how the specific choice of {rj} 
will influence Anderson localization effects. In 3d this 
is rather a number theory problem which is beyond the 
scope of this paper. Here in order to qualitatively address 
this issue we have employed two different methods to 
generate a given set of {n}: 



Method 1: We set t\ — a/ A 2 , r 2 = a/ A and r 3 
i\ where A = 1.3247.... is the real root of the equation 
A 3 — A — 1 = 0, and a is a random number extracted 
from a box distribution with support (10,20). In this 
way Ti 3> 1. This is necessary condition in order to have 
pseudo-random features. A given a is only accepted if 
the resulting continuous fraction form of tx/'k, T2/n and 
ts/tt specified by a sequence of integers, say, [0,1,0,2,—] 
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is such that the maximum value of |<Zj| is smaller than 
12. This is enough to prevent the system to be close to 
a quantum resonance. We note that the smaller these 
integers at are the more irrational are ti/tt. 

Method 2: We set t\ — a/ a, t 2 = a and T3 = act 
where a rs 1.618 is the golden mean. Then we follow the 
Method 1. 

In this section we carry out a finite size scaling analysis 
in order to determine the impact of this type of pseudo 
randomness on Anderson localization effects. We will 
show that the system undergoes a metal insulator tran- 
sition at k = k c = 2.27 ± 0.07 if Method 1 is used and 
k = k c = 2.35 ± 0.1 if Method 2 is used. By contrast if 
the phase of the kinetic term of the evolution matrix is 
completely random (see last section) k c = 2.40 ±0.05. It 
thus clear that the pseudo random nature of the {t^} have 
a sizable but small effect on the localization transition. 
We shall see in the last section that for certain choices of 
{Tj} these differences can be much larger and likely lead 
to k c quite different from the one corresponding to the 
random case. 

In order to proceed we evaluate the quantum evolution 
operator U. After a period T = 1 (we remember that 
h = 1), an initial state ipo evolves to 

V>(1) = Uifo = e^V^Vo, (15) 

where p and q stand for the usual momentum and posi- 
tion operator. We then express the evolution operator U 
in a basis of the momentum eigenstates, 

(m\U\n) = L e -i(rml+r 2 ni+r 3 nl)/2 e **&™*) (16) 

1 

where (q|n) = { f^ /2 with q = (qi,q 2 ,q 3 ), n = 
(ni.fte.na), <Kl,m,n) = 27r(l + 0)(m-n)/iV-F(27r(l + 
0)/N), 1 = {h, h, h), h = -N/2, . . . , A/2 - 1. Here N 3 
specifies the dimension of the Hilbert space. We restrict 
ourselves to even N throughout this work. The param- 
eter 9 depends on the boundary conditions < 0, < 1. 
We set 8 = that corresponds to the periodic boundary 
conditions. 

The eigenvalues and eigenvectors of U can now be com- 
puted by using standard diagonalization techniques. The 
eigenvalues of the evolution matrix are the starting point 
to carry out the finite size scaling analysis. 

Finally we note that following the procedure of Ref. 
Q it is possible to map a given eigenstate of U with 
eigenphase e luJ onto a 3d AM, 

m 

where 

e n = tan(w/2 - nn 2 /4 - T 2 n 2 2 /A - T Z n\j '4) (18) 

and 

F(m - n) = -7-^ [ d 3 gtan(^(q)/2)e- icl ( m - n ) 
( 27I T J-« 

oc exp(— yl|m — n|) 



with A a decreasing function of the kicking strength. The 
form of the diagonal disorder and hopping is similar to 
the Id case. We thus expect that the quantum dynamics 
of Eq. (Ti"T|) is qualitatively similar to that of a 3d AM with 
uncorrelated disorder and hopping restricted to nearest 
neighbors. However deviations are expected since the 
potential Eq. (TT8"| have strong correlations (see Fig. [5]). 



A. Finite size scaling method 

Our first task is to find out the critical k c and the crit- 
ical exponent v associated with the divergence of the lo- 
calization length at the transition. In order to proceed we 
determine k c by using the finite size scaling method [20| . 
First we evaluate a certain spectral correlator for differ- 
ent sizes N and disorder strengths k. Then we locate the 
transition by finding the kicking strength k c such that 
the chosen scaling variable becomes size-independent. In 
our case we investigate the level spacing distribution P(s) 
[25I ] (probability of finding two neighboring eigenvalues at 
a distance s = (e^+i — ej)/A, with A being the local mean 
level spacing). In order to avoid any dependence on bin 
size, the scaling behavior of P(s) is examined through 
the following function of its variance [26] 

r)(N, W) = [var(s) — varwo]/[varp — var-yvD]- (19) 

In Eq. (fT^j) var(s) = (s 2 ) — (s) 2 , where (...) denotes 
spectral averaging over a single set of {t^}. varwD ~ 
0.286 and varp = 1 are the variances of Wigner-Dyson 
(metal) and Poisson (insulator) statistics [25[ , respec- 
tively. Hence 77 = 1(0) for an insulator (metal). Any 
other intermediate value of rj in the N — * 00 limit is an 
indication of a metal insulator transition. 

In Fig. 1 we plot the k dependence of rj for different 
system sizes. The critical disorder k = k c signaling the 
Anderson transition corresponds to the point for which 
7] is independent of N. For a smaller (larger) k, rj tends 
to the insulator (metal) prediction. For a precise deter- 
mination of the critical k c and the critical exponent v we 
look at the correlation length near k Cl 

m = &\k - k c \- v , (20) 

where £0 is a constant. The numerical values of k c and v 
are obtained by expressing rj(N, k) = f[N/^(k)] and then 
performing an expansion around the critical point 

n(N, k) =Vc + J2 C ^ k - kc) n N n ^ . (21) 

n 

In practice, we have truncated the series at n = 2 
and then we have performed a statistical analysis of the 
data with the Levenberg-Marquardt method for nonlin- 
ear least-squares models. The most likely fit is deter- 
mined by minimizing the x 2 statistics of the fitting func- 
tion HD. We found k c = 2.27±0.07 with v = 1.61 ±0.25 
for Method 1 and k c = 2.35 ±0.1 with v = 1.67 ±0.27 for 
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FIG. 2: The spectral rigidity A3 (A) and level spacing distri- 
bution P(s) for three different kicking values, above, below 
and close to the transition k c pa 2.4 in the 3d kicked rotor Eq. 
(|f 3|) with {li} the same as in Fig. 1(a). A transition from Pois- 
son to Wigner-Dyson statistics for a system with time reversal 
invariance (GOE) is clearly observed as the kicking strength 
crosses the critical value k ~ k c from below. For k ~ k c 
the level statistics have all the signatures of a metal-insulator 
transition such as a) an asymptotically linear A3 (A), b) level 
repulsion for s — + as in a metal (lim^o P( s ) = 0) and c) 
exponential decay of P(s) for s > 1 as in an insulator. In 
agreement with the OPT, level statistics are scale invariant 
for k « k c - 



Method 2. The numerical value of v is higher but quite 
close to the one found in the 3d AM at the transition 
v w 1.5. We stress that these results were obtained with 
only the data of Fig. 1, namely, an average over only two 
sets of {Tj} for each method. 

These findings are consistent with the heuristic picture 
that correlations make the potential less random and con- 
sequently weaken localization effects. We note that this 
conclusion (see next section) might be misleading due to 
the small number of utilized. Indeed in the final sec- 
tion we shall see that our results are typical, namely, for 
the majority of {r^} we shall get similar results. However 
there are some choices of irrational {r^} that can lead to 
larger deviations from the results presented in this sec- 
tion. 



II. LEVEL STATISTICS 

As was mentioned previously the analysis of level 
statistics is a powerful tool to investigate the metal in- 
sulator transition. At the transition level statistics are 
intermediate between that of a metal and that of a in- 
sulator. For instance, a) P(s) ~ s for s <C 1, a typi- 
cal feature of a disordered metal, is still observed at the 
transition; b) A 3 (A) = xA/15 for A > 1 with x < I 
[l7l l2ll I22I l23l U% . This is similar to the result for an 
insulator A3 (A) = A/15; c) likewise P(s) oc e~ As with 
s ^> 1 and A > 1. For an insulator P(s) = e~ s . 

We explore to what extent these features are also ob- 
served in our 3d kicked rotor. As is shown in Fig. [21 
qualitatively the level statistics of the eigenphases of the 
evolution matrix at k fts k c are strikingly similar to those 
of a 3d AM at the transition. The spectral rigidity A3 (A) 
is linear, there is level repulsion for s<l. the spectrum is 
scale invariant and P(s) decays exponentially for s>l. 
The slope of A 3 (A), x ~ 0-29 in our case, is larger than 
the one corresponding to a random kinetic term (see next 
section). This suggests that the quasi-random nature of 
our model only modifies quantitatively, non-qualitatively, 
the transition. As was expected, level statistics in the 
limits k ^> (<C)fc c are described by Wigner-Dyson (Pois- 
son) statistics typical of a metal (insulator). We note that 
for systems with time reversal invariance the Wigner- 
Dyson statistics is usually defined by the spectral prop- 
erties of the Gaussian Orthogonal Ensemble of random 
matrices (GOE). 

III. QUANTUM DIFFUSION 

We now study quantum diffusion in the region k « k c . 
Our motivation is to find out whether the predictions of 
the scaling theory still hold in the 3d kicked rotor Eq. (fF3]) . 

Provided that the classical phase is fully chaotic we 
expect the classical motion to be diffusive (Ap 2 (t)) oc t. 
By contrast quantum dynamics depends strongly on k. In 
analogy with a prediction of the OPT for the 3d AM, we 
expect dynamical localization for k < k c . In the opposite 
limit, quantum effects are small, dynamical localization 
is avoided and diffusion does not stop though it might 
slow down as a consequence of interference effects. 

According to the OPT (see introduction) for k w k c , 
quantum diffusion must be anomalous {Ap 2 (t)) oc i 2 / 3 . 
In Fig. [3l we show the classical and quantum (Ap 2 (t)) 
for different values of the kicking strength. We find full 
agreement with the predictions of the OPT: dynamical 
localization is avoided for k > k c and, for k « k Cl diffu- 
sion is anomalous with {Ap 2 (t)) oc t@ , (3 w 0.66. We note 
that the determination of k c using the finite size scaling 
method is fully consistent with the results from quantum 
diffusion. 

The classical (Ap 2 (t)) was computed by averaging over 
an ensemble of 10 6 initial (t = 0) states uniformly dis- 
tributed in — 7r < 3j(Q) < it with Pi(0) = Q;i — 1,2,3. 
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FIG. 3: Classical (a) and quantum (b) diffusion, {Ap(t) 2 } = 
(\p(t) — p(0)] 2 }, for three different kicking values, above, below 
and close to the transition k c ~ 2.4 in the 3d kicked rotor 
Eq. (|13p with {r^} as in Fig. 1(a). Classical diffusion is normal 
for the three kicking values however interference effects slows 
down and eventually arrests quantum diffusion for k < k c . In 
agreement with OPT {p 2 ) oc t 2 ^' < at k ~ k c . Thin solid lines 
in (b) are the best fits in the window 10 < t < 1000. 



In the quantum case we average over the following four 
initial conditions, |n) = |0, 0, 0), |0, 0, 1), |0, 1, 0) and 
|1,0,0). 



FIG. 4: Upper: r\ versus the kicking strength for different 
system sizes N for a 3d kicked rotor in which the phase 
of the kinetic term of the evolution matrix is random (see 
text). A transition is again observed but for a slightly big- 
ger k c = 2.4 ± 0.05. Lower: r\ for a randomized 3d kicked 
rotor (random phase plot) and for the deterministic Eq. (|13[) 
one with {n} chosen from the two methods mentioned in the 
text. In both cases N = 16. Small but sizable differences are 
clearly observed specially in the insulator region. For these 
\t{\ correlations in the deterministic 3d kicked rotor seem to 
slightly weaken localization effects. In order to a meaning- 
ful comparison in the deterministic case we carry out a full 
ensemble average over 1280 sets of {n}. 



IV. EFFECT OF CORRELATIONS: 
COMPARISON WITH A 3D KICKED ROTOR 
WITH RANDOM KINETIC TERM 

In this section we compare previous results with those 
coming from a 3d kicked rotor similar to Eq. (|13p but 

with the kinetic term of the evolution matrix completely 

'i 2 _i_ 2 j- 2 ^ 
random, namely the factor e HTlTl i +T2Tl2+r3 "3-' is replaced 

by e*^ where <j> is a random number in [0, 2w). 

Our motivation is to clarify to what extent the deter- 
ministic but pseudo-random character of the motion of 
the 3d kicked rotor Eq. (fT3")) weakens Anderson localiza- 
tion effects. This is relevant for applications as in realistic 
situations it is hard to produce a truly random potential 



free from correlations. In essence we repeat the calcu- 
lations of the previous section but for the randomized 
version of the 3d kicked rotor: a) we carry out a finite 
size scaling analysis to determine k c and v (see Fig. 0}, 
b) we compute the spectral rigidity and level spacing dis- 
tribution in the metallic, critical and insulator region, c) 
we study quantum diffusion for k k c . 
The main conclusions are as follows: 

a) A transition is also observed in this model but at a 
slightly greater k ~ k c = 2.4 than in the deterministic 
case (see fig HI) ■ 

b) The slope of the number variance \ ~ 0.22 is smaller 
than in the deterministic case. This is consistent with 
weaker localization effects in the deterministic case. 

c) v w 1.6 ± 0.2 is roughly the same in both cases. This 
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FIG. 5: Distribution of jy for a 3d kicked rotor Eq. (|13[) with (a) the kinetic term replaced by a random phase, (b) {n} given 
by Method 1, (d) {n} given by Method 2, and (c) the same for Method 2 but where the integer sequence [01,02,...] for the 
continuous fraction form of n/n satisfies \ai\ < 9 instead. In all cases 1280 sets of {n} are considered and N = 16. 



is expected as the numerical value of critical exponents 
should not depend on the microscopic details of the sys- 
tem. 

d) Quantum diffusion is anomalous with (p 2 (t)) cx i 2 / 3 
and virtually indistinguishable from both the determin- 
istic and the 3d AM. This together with the scale invari- 
ance of the spectral correlations, the level repulsion, the 
asymptotic linear behavior of the spectral rigidity and 
the numerical value of critical exponent v seem to be the 
more robust features of the transition. 

A remark is in order. The results for the random case 
were obtained by performing an ensemble average over 
1280 different realizations. By contrast in the determin- 
istic case we did not want to perform an artificial ensem- 
ble average so the results presented corresponds to only a 
handful of different A natural question to ask is to 

what extent the numerical values of k c will be different 
if other {Tj} would have been used. For the {t^} used 
here we have obtained a k c which is slightly smaller than 
in the random case. By contrast in a recent paper [28[ 
we found k c w 3.3 and ?7 C « 0.06 for t\ = 1, T2 = 1/A 
and T3 = 1/A 2 with lambda the root of x 3 — x — 1 = 0. 
This value of k c is sizable larger than the one k c « 2.4 
found in the random case. This suggests that for that 



choice of {t;} the effect of correlations enhances local- 
ization effects. Several questions then arise: what is the 
probability that a given set of periods {t{\ leads to an en- 
hancement/supression of localization?, on average, what 
is the magnitude of the typical deviations from the ran- 
dom case? 

Unfortunately a detailed answer to these questions 
faces technical difficulties as the finite size scaling method 
is quite expensive numerically. Moreover, on the more 
fundamental level, this is rather a number theory prob- 
lem that, as far as we know, have not yet been tackled in 
the mathematical literature. 

In order to give a qualitatively answer to these ques- 
tions we compute the full distribution of ?y, its average 
(77) and standard deviation a v where (...) in this case 
stands for spectral and ensemble average (see Fig. 5). 
In the deterministic case we carry the ensemble average 
over different sets of {n}: from Method 1, Method 2 and 
Method 2 but with |a»| < 9. In each case 1280 differ- 
ent choices of {r^} are considered. In the random case 
the average is carried out also over 1280 different ran- 
dom realizations as it is explained above. Several con- 
clusions can be extracted from this numerical calculation 
(see Fig. 5): a) The distribution in the random case is 
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much narrower than in the deterministic case and as a 
consequence the standard deviation is much smaller. It 
is thus expected that specific choices of {r^} could lead to 
sizable deviations between the random and deterministic 
case, b) The distribution of 77 for the different methods 
of obtaining {7^} are qualitatively similar. That suggests 
that provided that is sufficiently irrational (\di\ < 12) the 
distribution of 77 is not very sensitive to the way in which 
{rj} is generated, c) The distribution of 77 is clearly asym- 
metric. A longer tail is observed in the region of greater rj 
which indicates that for a sizable percentage of {t,} local- 
ization effects are considerably suppressed. At the same 
time the maximum of the distribution is slightly shifted 
to the left with respect of the random case. Therefore 
for a typical {n}, k c will be somehow smaller than in the 
random case and consequently localization effects will be 
mildly suppressed due to the non random character of 
the motion. 



V. CONCLUSIONS 

We have studied the quantum dynamics of a 3d kicked 
rotor with a smooth potential. A careful finite size scaling 
analysis has shown that this system undergoes a metal- 
insulator transition for k « k c . The value of k c barely 



depends on the numerical choice of {tj} provided it is 
sufficiently irrational. However from the analysis of the 
distribution of 77 it is likely that for particular choices of 
irrational {r^} stronger deviations are observed. 

In the critical region k « k c we have investigated in 
detail quantum diffusion, level statistics and critical ex- 
ponents. It has been shown that quantum diffusion, in 
agreement with the prediction of the OPT, is quantita- 
tively similar to that of a 3d AM at the transition. Typ- 
ical signatures of a metal-insulator transition such as a 
linear spectral rigidity or level repulsion have also been 
observed in our 3d kicked rotor for both deterministic 
and random kinetic terms. 
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